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Fermion mass generation in the strongly coupled U(l) 
lattice gauge theory with fermion and scalar fields of equal 
charge is investigated by means of numerical simulation with 
dynamical fermions. Chiral symmetry of this model is broken 
by the gauge interaction and restored by the light scalar. We 
present evidence for the existence of a particular, tricritical 
point of the corresponding phase boundary where the contin- 
uum limit might possibly be constructed. It is of interest as a 
model for dynamical symmetry breaking and mass generation 
due to a strong gauge interaction. In addition to the massive 
and unconfined fermion F and Goldstone boson tt, a gauge 
ball of mass ms — l/2mF and some other states are found. 
Tricritical exponents appear to be non-classical. 



I. INTRODUCTION 

Attempts to construct a theory with dynamical break- 
ing of global chiral symmetries in four dimensions, which 
could explain or replace the Higgs- Yukawa mechanism 
of particle mass generation, usually lead to the intro- 
duction of a new strong gauge interaction beyond the 
standard model and its standard extensions. For ex- 
ample the heavy top quark and the idea of top conden- 
sate 01 inspired the strongly coupled topcolor and simi- 
lar gauge models [|| (for a recent overview see e.g. [H ) . 
Among the requirements such a theory should satisfy, the 
most general ones are the following two: First, because 
gauge theories tend to confine charges in a regime where 
they break chiral symmetries dynamically, the physical 
states, in particular fermions, must be composite singlets 
of the new gauge symmetry. Second, as a strong coupling 
regime is encountered, the models should be nonpertur- 
batively renormalizable in order to be physically sensible 
in a sufficiently large interval of scales. 

Even in very simplified models, this are too difficult 
dynamical problems to get reliably under control by an- 
alytic means only. Therefore, a numerical investigation 
on the lattice of some prototypes of field theories with 
the above properties may be instructive. In such an ap- 
proach, the presumably chiral character of the new gauge 
interaction and numerous phenomenological aspects have 
to be left out of consideration. 

A promising candidate for such a prototype field the- 
ory on the lattice, the xU(f) model, has been described 
in Ref. Here the four-dimensional vector-like U(l) 



gauge theory contains the staggered fermion field x ^-nd 
the scalar field 0, both of unit charge. A Yukawa coupling 
between these matter fields is prohibited by the gauge 
symmetry. The global U(l) chiral symmetry, present 
when the bare mass mo of the fermion field x vanishes, 
is broken dynamically at strong gauge coupling g by the 
gauge interaction, similar to QCD or strongly coupled 
lattice QED. Whereas both x £^nd (f) constituents are 
confined, the massive physical fermion F = (f^x with 
shielded charge appears. 

The scalar suppresses the symmetry breaking when it 
gets lighter and induces a phase transition to the chi- 
ral symmetric phase [^,^. At this transition, for large 
enough gauge coupling, the mass of the physical fermion 
in lattice units scales, amp — > 0, and the lattice cutoff 
1/a thus can be removed for mp fixed in physical units. 
If the theory were renormalizable, a continuum theory 
with massive fermion as well as a massless Goldstone 
boson ("pion" tt) would be obtained. When the global 
U(l) chiral symmetry, modelling the SU(2) symmetry of 
the standard model, is gauged, this tt boson is "eaten" 
by the corresponding massive gauge boson. This is what 
is achieved in standard approaches by the Higgs- Yukawa 
mechanism. 

In this paper we address the question of renormaliz- 
ability of the xU<t> model at the line of chiral phase tran- 
sitions induced by the scalar field. We have no definite 
answer, but our extensive numerical study of the model 
in the relevant region of the three-dimensional parame- 
ter space (see Fig. |^) with dynamical fermions provided 
several encouraging results: 

1. Our previous studies have indicated that on 
the nearly whole chiral phase transition line, start- 
ing at the strong coupling limit (3 ~ g^^ = 0, the 
model behaves like the Nambu-Jona-Lasinio (NJL) 
model, belonging presumably to the same univer- 
sality class. We now present strong evidence that 
aX (3 — Pe — 0.62 the line contains a special point, 
the tricritical point E, where for theoretical rea- 
sons the scaling behavior is different from the rest 
of the line. It is governed by another fixed point. In 
difference to the NJL model the gauge field is not 
auxiliary but plays an important dynamical role at 
the point E. 
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FIG. 1. Schematic phase diagram of the xU4' model. Three critical lines, NE, EooE and E_ooE meet at the tricritical point 
E. The line NE is a part of the boundary of the Nambu phase (dark shaded region) at mo = with spontaneously broken 
chiral symmetry, amp 7^ 0. This phase is a sheet of first order phase transitions at which the chiral condensate changes sign. 
The lines EooE and E_ooE are critical boundaries of the "wings" of first order Higgs phase transitions. The light shaded region 
at mo = corresponds to vanishing fermion mass, amp ~ 0. The vertical sheets containing the points T and C separate the 
confinement and Coulomb phases. The line ETS is a line of tripple points. 



2. Using advanced methods of the finite size scaling 
analysis we estimate several tricritical exponents 
determining the scaling behavior at the point E and 
find nonclassical values. This indicates that due to 
the strong gauge interaction this point differs from 
the standard expectations for tricritical points in 
four dimensions |9|. 

3. In the vicinity of the point E in the broken phase, 
not only the fermion mass am p , but also the masses 
of several bosons (neutral states composed of scalar 
and gauge fields) and mesons (xx states) scale, i. e. 
in lattice units they approach zero with constant 
ratios. This suggests a rich spectrum if the contin- 
uum limit of the model is approached at the point 
E. In particular, a gauge ball of mass ms — 1/2 mp 
is observed. 

4. The composite Goldstone state tt with properties 
required by chiral symmetry breaking is present. 

5. We determine the effective Yukawa coupling j/r be- 
tween the composite F and tt states in the vicinity 
of the critical line and find that lines of constant 
2/R tend to approach the point E [Q. We cannot 



yet say whether some of them end at this point, 
which would imply a nontrivial continuum limit. 
However, this approaching means that the coupling 
decreases only slowly with an increasing cutoff 1/a 
on paths towards E, thus increasing the chances for 
rcnormalizability. 

We have not been able to achieve at least qualitative 
results in two issues of major interest: A heavy scalar a- 
meson, which would correspond to a composite Higgs bo- 
son, is seen, but its mass in lattice units does not yet scale 
on the lattices of sizes we could afford and is strongly 
dependent on the bare fermion mass amQ. We cannot 
say anything about its value in the continuum limit at 
E. Also the pion decay constant /tt does not scale, i.e. 
f-n/TTLp seems to increase with decreasing distance from 
E. Its current value (at arriQ ~ 0.4) is about 1/3. The 
present data are consistent both with the possibility that 
/tt diverges in physical units, which would indicate triv- 
iality and that the absence of scaling is due to too 
small lattices. 

Concerning the possible triviality, we point out again 
[Q that the xU4> model would be a valuable model even 
if the cutoff cannot be removed completely without loos- 
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ing the interaction, provided the cutoff dependence of the 
renormahzed couphngs is sufhciently weak, e.g. logarith- 
mic as in the standard model. A dynamical approach 
to the Higgs- Yukawa mechanism does not necessarily re- 
quire a nontrivial fixed point. The Higgs- Yukawa sector, 
whose validity is restricted due to the triviality by a cer- 
tain upper energy bound, can be replaced by a theory 
with a higher upper bound. 

Nevertheless, it is possible that the xUtf) model in the 
continuum limit taken at the tricritical point E defines 
an interacting theory. The pursuit of this question re- 
quires a better understanding of tricritical points in four 
dimensions, as the available experience with such points 
is restricted to lower dimensions 0. Further obstacles 
are the necessity to tune two couplings and the need to 
extrapolate to the chiral limit, mo 0. Finally, more 
insight is needed into strongly coupled and not asymp- 
totically free nontrivial four-dimensional gauge theories, 
whose existence has been recently suggested by numerical 
investigation of pure U(l) gauge theory fl^ ]. 

We remark that the properties of the xt^0 model in 
lower dimensions are much better accessible. In two di- 
mensions, the numerical evidence strongly suggests that 
the continuum limit of the model is equivalent to the 
two-dimensional chiral Gross-Neveu model, and is thus 
renormalizable and asymptotically free jl^ . First results 
in three dimensions |Q suggest that the xUfj) model be- 
longs to the universality class of the three-dimensional 
chiral Gross-Neveu model, which has a non-Gaussian 
fixed point. In both cases the continuum limit is ob- 
tained on a whole critical line of chiral phase transitions 
emerging from the corresponding Gross-Neveu model ob- 
tained in the limit of infinite gauge coupling, without any 
use of possible tricritical points. In this sense the situ- 
ation in four dimensions is unique, and the experience 
from lower dimensions is not applicable. 

If the tricritical point E in the four-dimensional xC^0 
model defines a renormalizable continuum theory, simi- 
lar property might be expected in analogous models with 
other gauge symmetry groups. For example, an SU(2) 
gauge model with scalar and staggered fermion field in 
the fundamental representation of SU(2) is known to 
have at strong coupling a phase structure very similar 
to the xUff) model with the U(l) gauge field [^j6). There- 
fore we expect that the model we are studying is generic 
for a whole class of strongly coupled gauge models with 
fermions and scalars in the fundamental representation. 

After describing the xUfj} model in the next section, we 
present our results as follows: Some preparatory studies 
of the model in the limit of infinite bare fermion mass 
are presented in Section III. In the following section we 
demonstrate the existence of the tricritical point. In Sec- 
tion V the critical and tricritical exponents are estimated 
by finite size scaling studies. Spectrum in the continuum 
limit taken at the point E is discussed in Section VI. Then 
we summarize our results and conclude. In the Appendix 



we give a detailed definition of the meson propagators 
and effective Yukawa coupling we have calculated. 

Preliminary results of this work have been presented 
in Refs. ||l^,|l6|. An account of our results for the ef- 
fective Yukawa coupling between tt and F is given in a 
separate paper [Q. An investigation of the xUcj) model 
in the quenched approximation, with particular emphasis 
on the role of magnetic monopoles, has been performed 
in Ref. [0. A detailed presentation of the xU4' project 
in two, three, and four dimensions can be found in [ p^ . 
An investigation of similar models in continuum has been 
performed by Kondo . 

II. THE xU(l) MODEL 
A. Action and phase diagram 

The four-dimensional lattice x^ip model is defined by 
the action 

S^S^ + Su + S^, (2.1) 

X /.i— 1 

+amoJ2xxXa:, (2.2) 

X 

^c/ = -/3^cos(ep), (2.3) 
p 

4 

5^ = -K^^(0t[/^_^0^^^ + /j.c.). (2.4) 

X /I— 1 

Here 9p <E [0, 27r) is the plaquette angle, i.e. the argu- 
ment of the product of U(l) link variables Ux,^ along a 
plaquette P. Taking Op = a^gFf^,^, where a is the lattice 
spacing, and (3 = one obtains for weak coupling g 

the usual continuum gauge field action Sjj = j J d'^xF^^. 
The staggered fermion field x has (real) bare mass arriQ 
in lattice units and corresponds to four fermion species 
in the continuum limit. The scalar field (j) is of fixed 
modulus, |(/>| = 1. 

The model has U(l) global chiral symmetry in the limit 
Too 0, where toq is the bare fermion mass in physical 
units, to be defined while constructing the continuum 
limit. This is to be distinguished from the limit arriQ — > 
0, allowing explicit chiral symmetry breaking, toq ^ 0, 
when a ^ 0. Because of this fine difference between toq 
and aTOo, important in various possible continuum limits, 
we keep trace of a throughout the paper. 

The schematic phase diagram is shown in Fig. |^. We 
recognize several limit cases of the xU4> model as models 
interesting by themselves: 

1. At K = and atoq — oo, the pure U(l) gauge theory 
with the Wilson action (|2.3|) and phase transition 
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between the confinement and Coulomb phases. Its 
continuum limit in an extended coupling parameter 
space may be determined by a non-Gaussian fixed 
point [H). 

At K = and amo finite, t he g auge theory with 
fermions, i.e. compact QED (2.2) and (2.3), whose 
phase transition is currently under investigation 



mo approaches then infinity and fermions decouple. The 
remaining U(l) Higgs model is equivalent to the trivial 
theory at the critical endpoint of the Higgs phase tran- 
sitions ]2^ , p^ . This is confirmed by some of our results 
presented below. 



B. Observables and numerical simulations 



3. At /? = 0, i.e. the gauge field being auxiliary, the 
Nambu-Jona-Lasinio (NJL) model, obtained by in- 
tegrating out the bosonic fields The triviality 
of this model has been recently confirmed in large 
scale simulations . 

4. At atriQ = oo and k arbitrary, the compact scalar 
QED or U(l) Higgs model and (|^). Its con- 
tinuum limit at strong gauge coupling is Gaussian 
[§,11 ■ 

At strong coupling, /3 < 1, the model has three sheets 
of first order phase transitions: the two "wings" at finite 
amo, and the sheet at amo = 0, separating the regions 
with nonzero chiral condensate of opposite sign. These 
three sheets have critical boundary lines Ej-oqE and NE, 
respectively. As we shall discuss below, we have veri- 
fied with solid numerical accuracy that these 2"^^ order 
phase transition lines do indeed intersect at one point, 
the tricritical point E. We arc not aware of a convincing 
theoretical argument why this should be so. 

Of most interest is the Nambu phase at mo = 0, at 
small P and k. Because of confinement, there is no (f>- 
boson, i. e. charged scalar, neither fundamental charged 
X-fermion in the spectrum. The chiral symmetry is dy- 
namically broken, which leads to the presence of the neu- 
tral composite physical fermion F = (p'^x with the mass 
amp > 0. It scales, amp \ 0, when the NE line is 
approached. 

Further states include the "mesons" , i. e. the fermion- 
antifcrmion bound states: the Goldstone boson tt with 
am^ c>c y/amo, the scalar a, and the vector p. "Bosons" , 
observed in the scalar-antiscalar or gauge ball channels 
are present too. It is in particular the neutral scalar 
boson S. In the vicinity of the Ej-ooE lines the same 
scalar appears both in the (f>^ — (f) and gauge-ball channels 
which strongly mix. In the Nambu phase it is natural to 
interpret S* as a gauge ball, as this interpretation holds 
also deep in the Nambu phase, when the charged scalar 
(j) is heavy. 

The mass ams of the 5-boson vanishes on the lines 
EicxjE, whereas amp vanishes on the line NE. Both of 
them vanish at the tricritical point E. As their ratio is 
finite, the continuum limit obtained when approaching 
the point E contains both states and is thus different 
from the rest of the NE and E±ooE lines. 

The critical lines Ej-ooE provide another approach to 
the continuum. For a nonvanishing amo the bare mass 



For the investigation of the tricritical point we use the 
following observables: 

To localize the Higgs phase transition we use the nor- 
malized plaquette and link energy are defined as 



a;, iJ.<y 



(2.5) 
(2.6) 



where \=L^T is the lattice volume. Following [^,0 we 
use the perpendicular and parallel components of these 



energies. 



E±_ = Ei^ cos 9 + Ep sin 9 , 
i?l I = i?L sin 9 — Ep cos 9 , 



(2.7) 
(2.8) 



where 6 is the slope of the Higgs phase transition line at 
the endpoint in the plane (/3,4/3k). 

For the localization of the chiral phase transition we 
measure the chiral condensate 



(xx> = (TrA/-i) 



(2.9) 



via a stochastic estimator, where M is the fermion ma- 
trix. 

To calculate the mass of the physical fermion we con- 
sider the gauge invariant fermionic field Fx = (j^Xxx- The 
mass amp is measured by fitting its propagator in mo- 
mentum space 0. The results for the measurement in 
configuration space are consistent. 

The fermion-antifermion composite states are called 
"mesons" . The corresponding operators and other de- 
tails are given in [^|JT^. We tried to include also the 
annihilation part, but failed to obtain sufficient statis- 
tics. 

To improve the signal, we also measure the meson 
propagators with smeared sources. This required the 
adaption of the routines, used with Wilson fermions, to 
the case of staggered fermions. It is described in Ap- 
pendix A 1 . With these smeared sources we have been 



able to fit the meson propagators by a one particle con- 
tribution at time distances larger than zero. But the 
same masses could be obtained if the unsmeared prop- 
agators were fitted with the inclusion of excited states. 
The smeared propagators reduce the errors, however, and 
in this work we mostly show results obtained by this 
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method. Further details of the fitting procedure can be 
found in [T^ . 

From the propagator of the n meson also the pion de- 
cay constant a/tt can be calculated pj], 



anio 



Here am,r and are the mass and the wave function 
renormalization constant of the tt meson. We checked, 
that afrr fulfills with excellent precision the current alge- 
bra relation 



(2.11) 



This is so even very close to the phase transition, though 
both o/tt and (xx) show rather strong finite size effects 
there. 

For the investigation of the chiral phase transition we 
also calculate the susceptibility ratio i?7r, which is de- 
fined as the logarithmic derivative of the chiral conden- 
sate l25|,ra, 



Rtt — 



91n(xx) 



d In arriQ 



ama d{xx) 



/3,K 



(XX> damo 



(2.12) 



/3,K 



We measure it as the ratio of zero momentum meson 
propagators 



a(p^o) 
a(p = o) 



(2.13) 



including the annihilation part of the propagator. This is 
done by means of a stochastic estimator and is described 
in detail in . 

As explained in we expect that close to a critical 
point the data could be described by means of the scaling 
law 



where t is the distance from the critical point (reduced 
coupling), A = p+j the critical exponent and G a scaling 
function. At the critical point, 



Rtt{0, amo) = g{oo) 



1 

J' 



(2.15) 



as ca n be seen by inserting (xx) oc {amo)^^^ into equa- 
tion (2.12). At the critical point, i?^ should be indepen- 
dent of amo for sufficiently small amo (scaling region). 
In the broken phase, i?^ vanishes in the chiral limit, as 
can been seen easily from the definition. In the symmet- 
ric phase, the a and tt channels are degenerate, so that 
in the chiral limit R^^ = 1. For small fixed t a charac- 
teristic behavior is expected, if one varies amo: Because 
A > 1, close to the critical point the curves for R-^ start 
for amo = from and 1 , respectively, and for increasing 



amo approach the horizontal line 1/5. This will happen 
the faster the smaller |t| is (compare fig. [T^ ). 

Further we consider the scalar and vector bosons, 
whose operators are defined as 



(2.10) 0(5)(,) 



-l5:Re|^4,t/(,^ 

X (. 4=1 



t),i(Px+i,t ( , 



(2.16) 



^r^W ^ I3 El°^ {4,tUis,tu<l^s+r,t} , (2.17) 
i = 1,2,3. 

The masses ams and amy of the scalar and vector bosons 
are calculated from the corresponding correlation func- 
tions in configuration spaceQ. 

In the same way we also measure the gauge invariant 
combinations of the gauge fields, which we call gauge 
balls, in analogy to QCD glue balls. We define two oper- 
ators with the quantum numbers 0"*""*" and 1"' : 



3 i-l 



o(°^^)(^)=i,5:ReK:^c/(,,, 

X \^i=2 j=l 

X 

i = 1,2,3. 



O 



(2.18) 
(2.19) 



The masses amQ of the gauge-balls are calculated in anal- 
ogy to the boson masses by means of the propagators in 
the configuration space. 

We have observed mixing of the S boson and the O^"*" 
gauge ball by means of the two-point function 



(2.20) 



We have also measured the effective Yukawa coupling 
j/R between the neutral fermion F and the tt meson. This 
is done in analogy to J2^, and the used operators are 



(2.14) described in appendix A 2. A detailed discussion of our 



results is given in pC| ] and summarized in the conclusion 
of the present paper. 



III. LIMIT OF INFINITE BARE FERMION MASS 

A. U(l) Higgs model and chiral phase transition in 
the quenched approximation 

For amo — oo, the xC^0- model reduc es t o the U(l ) 
Higgs model with |(/)| = 1 on the lattice, ( |2.3| ) and (|]j). 



^To reduce the statistical fluctuations in the determination 
of ams, calculating the propagator we subtract the momen- 
tum zero propagator before the determination of the error 
(average over the propagator). 
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This model has been investigated in the eighties (for a 
review see e.g. |^) and with modern methods in p2| , p3t . 
Its phase diagram is represented by the front face of Fig. 
|l|. It has the Coulomb phase at small n and large (3, the 
rest being the confinement-Higgs phase. The line of Higgs 
phase transitions EooSoo is first order except the points 
Eoo and Soo • The continuum limit at the critical endpoint 
Eqo corresponds most probably to a trivial scalar field 
theory ||,||. 

When dynamical fermions with amo > are included, 
the phase diagram remains roughly the same, except that 
the confinement- Coulomb phase transition and the end- 
point Eamo shift to smaller (3. The endpoints then form 
the critical line EqoE. It is natural to expect that this 
line, except the tricritical point E, remains in the same 
universality class as the point Eoq. Our results confirm 
this expectation. 

When quenched fermions with small amg are included 
into the Higgs model, a line of chiral phase transitions 
appears in the otherwise unchanged phase diagram of 
the Higgs model. It was realized already in the first in- 
vestigations of the Higgs model with fermion, that full 
and quenched models have a very similar phase diagram 
. This includes the observation that the chiral phase 
transition line runs within numerical accuracy into the 
critical endpoint of the Higgs phase transition line. The 
phase diagram of the quenched model looks thus similar 
to the arriQ = plane of Fig. 0. 

This similarity suggests that it might be instructive 
to study the xC/0-model in the quenched approximation. 
In Ref. ijl^ a quenched investigation of the interplay of 
chiral phase transition and the monopole percolation was 
performed. It seems that there might be an interplay of 
both transitions at an intermediate (3 on the NqqEoo line, 
possibly with nontrivial exponents. Around the points 
Noo and Eoo the chiral and the percolation transitions 
appear to be separated, however. 



B. Scaling behavior at the endpoint Eoo 

We begin by investigating the endpoint Eoo ■ We want 
to gain experience and check the reliability of the deter- 
mination of critical exponents by means of Fisher zeros. 
We later apply this method at finite amo for the scaling 
investigation along the EooE line and compare the results 
with those at Eoo. 

The scaling behavior at the endpoint of the Higgs 
phase transition line was determined in p^p^ along the 
first order Higgs transition line. It was found that the 
endpoint is described by mean field exponents. We inves- 
tigate the scaling behavior approaching i?oo in different 
directions. For this purpose it is useful to introduce the 
following reduced couplings (Fig. ^): 

t : parallel to the phase transition line of 1*^* order 
h : perpendicular to the phase transition line. 




FIG. 2. Definition of the reduced couplings t and h and of 
the angle 6. The point Eoo is the endpoint of the Higgs phase 
transition line. The scale of the k axis is 4/3 times larger than 
that of the /3 axis. 



Here perpendicular is understood in the same sense as in 
eqs.(|2!7|) and (|U), so that 



t 
h 



'{P — Pc) COS 6* -I- I (k — Kc) sin ( 



{(3 — f3c) sinO -I- |(k — Kc) cos( 
and therefore 

S = ~QV{tE\\+hE^) + SE^ ■ 



(3.1) 
(3.2) 



(3.3) 



The letters t and h have been chosen in analogy to tem- 
perature and external field in magnetic systems. 

We introduce critical exponents v and v of the corre- 
lation length for both directions, 



and 



\h=Q ; 



\t = ; 



(3.4) 



(3.5) 



^ being the correlation length diverging at Eqo. 

To understand the relation between v and v we assume 
the equation of state 



(3.6) 



with the scaling function F and A = /3 -I- 7. Introducing 
F{x) = x'^ F{x^) it can be rewritten as 



e = \h\-''/^F 



\h\ 



l/A 



1^1 



(3.7) 



Assuming F{po) < 00 this means D — v/A. 

The scaling behavior (3.5) is expected in the general 
direction t = c ■ h, because A > 1, and therefore 



\h\ 



l/A 



\t\ 



00 



(3.8) 



for h ^ and, accordingly, t ^ 0. This makes clear that 
it is not important to choose h perpendicular to t. Only 
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the t-direction {h = 0) is special as it is tangential to the 
phase transition line and thus described by the scaling 
law (U). 

The mean field values of the exponents [3 — 1/2, 7=1, 
and 1/ = 1/2 correspond to i> = 1/3. 

To determine the critical exponent of the correlation 
length we measure the scaling behavior of the edge sin- 
gularity in the complex coupling plane (Fisher zero) . 
From scaling arguments for the free energy we expect for 
the first zero zi: 



Imzi(L)|j^ 



A-L 



(3.9) 



As all directions which arc not tangential to the phase 
transition line are equivalent, we expect the same ex- 
ponent v also if we fix /J = p-^,^ or k = k-^^ . This was 
verified in []l5| . Fixing one of the couplings is particularly 
convenient for the necessary analytic extrapolations into 
the complex plane. That is done by means of the multi- 
histogram reweighting method [^0| . 

We present here the scaling investigation we did for 
P = 0.848 « /3eoo- Fig- ||a shows a nice scaling be- 
havior for all lattice sizes with the critical exponent 
i> = 0.3236(10). This value is very close the mean field 
exponent i> — 1/3. 

The small deviations outside the error bars are proba- 
bly due to logarithmic corrections, as they are expected 
at Gaussian fixed points. To verify this we follow the idea 
from ||3^ and factor out the leading power law L-^/"^. If 
the scaling behavior has the form 



lulZl{L)\^ 



A ■ (InL)- 



(3.10) 



we expect in the In In plot a strait line with the slope —p. 
The data shown in Fig. |^ are very well described by a 
straight line with p ~ 0.17. The i s sm aller than that 
with a fit by means of the equation ( |3.9[) . However, we 
have not investigated how far these results, and especially 
p depend on the precise knowledge of the critical point. 

A similar value for the exponent j> was also measured 
in the SU(2) Higgs model ||3^ , ^ , but then it was not re- 
alized that this value is actually compatible with a Gaus- 
sian fixed point. Now we can conclude that the U(l) and 
SU(2) Higgs models have a very similar scaling behavior 
at the endpoint. 

We have also determined the critical exponent iy by 
means of the finite size scaling of the specific heat and of 
some cumulants [T^ ] . These less precise methods confirm 
the results presented here. Also the shift exponent A has 
turned out to be compatible with l/i> for the investigated 
observables ITslJil. 



IV. EXISTENCE AND POSITION OF THE 
TRICRITICAL POINT E 



(a) 




fit: A'L""': 

0.32355(99) 
A= 0.760(11) 

);^„=0.B3I250 



(b) k 




p= 0.166(17) 
lnA= -0.3451(76) 
X/dof=0. 152078 



0.4 0.6 0.8 1.0 

ln{ln I,) 

FIG. 3. (a) Scaling behavior of the edge singularity for 
f3 = 0.848 ~ Pe^ in the U(l) Higgs model on L* lattices. 
The small deviations from the trivial scaling behavior with 
i> — 1/3 are resolved in (b). Here ln{L^Imzi) is shown as a 
function of In(lnL). The fit from Fig. (a) is shown dotted. 
The value for /3e^ has been taken from |23|. 



A. General properties of trier itical points 

To present our investigations of the tricritical point 
E, we first summarize the relevant general properties of 
tricritical points and define the exponents. In notation 
we follow Grifiiths 

In the vicinity of a tricritical point it is usual to choose 
the following orthogonal coordinate system (Fig. ^): 
A : tangential to the first order PT line in the symmetry 
plane, 

g : perpendicular to the PT line inside the symmetry 
plane, 

C : perpendicular to the symmetry plane. 
In the symmetry plane mg = 0, these definitions are anal- 
ogous to those in the Higgs model (sec. HIB ), A and g 
corresponding to t and h, respectively. 

In the phase diagram there are four special lines, which 
we denote following the chiral PT line NE (second 
order) in the symmetry plane (A > 0) is lambda line 
La, its continuation in the symmetry plane, on which 
three first order phase transition sheets meet (A < 0), 
is triple line L,-, and the two lines of endpoints outside 
the symmetry plane are wing critical lines L4- and L_. 
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FIG. 4. Schematic phase diagram of the first order PT 
planes and the critical lines (L±, L^) in the vicinity of the tri- 
critical point E. Lt is the line of triple points. The L lines cor- 
respond the lines in Fig. ^ as follows: La = NE, L± = E±ocE, 
Lt =ET. a local coordinate system is shown. 



The first order PT plane below the lines La and hr in 
the symmetry plane is denoted So, and the two wings of 
Higgs phase transitions are S+ and S_. Because of the 
±amo symmetry we use in the following only the index 
+. 

The most important exponents and the defining scal- 
ing behavior are summarized in table ^ Comparing a 
metamagnet to our model, the staggered magnetization 
corresponds to the chiral condensate and the magneti- 
zation to the energy E±. The unfortunate fact is that 
the symmetry plane, which is of major use in the study 
of metamagnets, corresponds to the chiral limit toq = 0, 
which is difficult to approach in numerical simulations 
with fermions coupled to a gauge field. 

The two sets of exponents with index t and u are de- 
fined in analogy to the exponents on the adjacent critical 
lines: The set with the index t (tricritical exponents) is 
defined in analogy to the exponents along the chiral PT 
line NE (L^ line). The set with the index u (subsidiary 
exponents) is defined in analogy to the exponents along 
the line of Higgs endpoints EooE (L_|_ line). In general, 
the exponents at the tricritical point are different from 
those at the adjacent lines. We denote the diverging cor- 
relation lengths and the exponents on the lines and 
L+ by the indices A and respectively. 

The tricritical point can be considered as a special 
point of both the L^ and the L+ critical lines. Both 
the correlation length diverging at the La line and the 
correlation length ^_|_ diverging at the L_|_ line are criti- 
cal there. Nevertheless, in general, and ^+ are to be 
distinguished also at the tricritical point. In our case 
^A = l/amp and ^+ = l/ams- 
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Aii = PuSu 
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TABLE I. Exponents at the tricritical point. In the second 
column the notation of |^ is given. The classical value is de- 
rived in three dimensions. Its applicability to four dimensions 
is discussed in the text. 



point E, ^A ^+ on all paths into the point E. This 
proportionality seems to hold also for all other observed 
correlation lengths (inverse masses) which diverge at the 
point E. Thus there seems to be only one scaling law and 
i^u = ^t- This is a generic property of tricritical points. It 
makes possible to use at the tricritical point the finite size 
scaling theory quite in analogy to the adjacent critical 
lines. 

In three dimensions it is usual that tricritical points 
have a large region of dominance. In analogy, near the 
tricritical point we expect to find at some distance from 
the second order PT lines already the scaling behavior 
described by the tricritical exponents. Such a crossover 
phenomenon was investigated for example in . A sim- 
ilar effect might be expected also on small lattices in the 
immediate vicinity to the PT lines. It is not excluded, 
however, that in four dimensions the tricritical points 
are much less dominant than in three-dimensional mod- 
els. To the best of our knowledge, tricritical points have 
not yet been investigated numerically in four dimensions. 

For the singular part of the free energy F usually the 
following scaling relations are assumed : 



i^sing(A,g,C) = |Ap-""^(±)(^ 
F,i„g(A,g,C) = |.9p-"'^P^ 



C 



lAI"^' |A|^^ 
A C 



(4.1) 
(4.2) 



Our results strongly indicate (sec. |VI A 21 ) , that at the 



isr'"^ \9\ 

For such systems several relations between the exponents 
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can be derived |34| : 

= 1 - , 

2 - a„ = (/)(2 - at) , 

l + 6t = [2- at)/ (it, 

At=(3t6t = {2~au)/^-6t/{l + St), 
Pt = {2-at)/{l + St). 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 



In our work we use in particular the last two of these 
relations. 

Only four exponents are independent. With the as- 
sumption that the hyperscaling relation 



at = 2 - di^t 



(4.10) 



holds, only three independent exponents remain. 

Unfortunately, theory of tricritical points in four di- 
mensions is developed only insufficiently. The Ginzburg 
criterion indicates that for all dimensions c? > 3, tricrit- 
ical point can be described by the classical exponents. 
This would correspond in four dimensions to a violation 
of the hyperscaling relation . 



B. Two diverging correlation lengths 

For the existence of the tricritical point E in the xU^i 
model the chiral and Higgs phase transition must meet 
at one point in the toq = plane. Since there is no 
theoretical understanding for the interplay between both 
transition, the existence of such a point has to be demon- 
strated numerically. 

To give a first impression. Fig. || shows the mass of the 
scalar boson ams in the vicinity of the tricritical point for 
aniQ = 0.01, 0.02, and 0.04. This mass h as be en obtained 
from the (j)'^ — (j) correlation function ( |2.16 ). It has a 
pronounced minimum (arrows in Fig. |^) for each aniQ. Its 
minimal value on the 6'^ 16 lattice is about 0.35, and 0.2 on 
the 8'^24 lattice. The significant decrease with increasing 
lattice size indicates that the mass in the infinite volume 
vanishes, and the correlation length = 1/ams thus 
diverges. The vanishing of ams at some (3, k for any finite 
aniQ corresponds to the EoqE line of cndpoints Eamo of 
the Higgs phase transitions. 

The fermion mass amp descends steeply (Fig. ^ at 
the position of the minimum of the boson mass for the 
same amo and volume. In the broken phase, the curves 
for different volumes slowly approach each other. In the 
symmetric phase, the values of amp achieve small finite 
values which should vanish in the chiral limit and infinite 
volume. 

Figures ^ and ^ show that changes of volume and amo 
result in a shift of the minimum in p. The same holds 
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FIG. 5. Boson mass ams as function of l3 for different amo 
and lattice sizes at k = 0.30 ~ ke. The arrows indicate the 
minima of ams. 
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FIG. 6. Fermion mass amp as function of (3 for different 
amo and lattice sizes at k = 0.30. The minima of the bo- 
son mass ams are marked with full symbols and arrows (cf. 
Fig. i). 



for K. So there is little hope to extrapolate the data at 
fixed P and n into infinite volume and chiral limit. As 
usual for tricritical points, a fine tuning of both couplings 
is required. In this work we assume that the limited 
precision of the position of the tricritical point we have 
achieved is outweighted by a sufficiently large domain of 
dominance of this point. 



C. Position of the tricritical point E 

To localize the point E we determine the positions of 
the endpoints J^amo for small amp and extrapolate them 
to amo = 0. It is difficult to control the uncertainties in 
every step of this procedure. 

At nonzero amo we proceed in analogy to p3[ . We 
determine the latent heat of E±_ in the S+ plane. Fig. |^ 
shows the latent heat AE± as a function of /3 for two 
amo and different lattice sizes. The values have been 
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FIG. 7. Square of the latent heat of the energy Ej_ for 
amo = 0.02 and amo = 0.04 for lattices of the size 4^ to 8* 
and 10'', resp.. The mean-field exponent (3+ = ^ corresponds 
to a strait line. We have desc ribed the data for each amo with 
one common /3+ in eq. (4.11) for the different lattice sizes. 
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FIG. 8. Pseudocritical couplings (3pc{L) as function of 
for different amo < 0.10. The fit is an extrapolation into 
infinite volume assuming the exponent value = 0.5. 



amo 



obtained by reweighting the data to the couphng where 
both maxima have the same height. Then the distance of 
these maxima in the histogram and the uncertainty have 
been estimated, as the data is not sufficient for a fitting 
procedure. 

We expect a scahng of AE±_ for fixed aniQ and fixed 
lattice size of the form 



AEj_ (X t/^^ (x{(3- f3pc{L, amo)) 



/5+ 



(4.11) 



The index + denotes the magnetic exponent on the EqoE 
Hne. In this procedure it is assumed that the dominant 
contribution due to the finite volume can be absorbed in 
a volume dependent pseudocritical coupling f3pc{L, amo). 
A better method would have been to extrapolate the la- 
tent heat first into the infinite volume and to investigate 
scaling afterwards. But for such an analysis the quality 
of our data is not sufficient. 

In the quadratic plot in Fig. ^ we expect a straight 
line for the classical value /3+ = ^ , which was observed 
for amo = oo p^ ]. Our data at small amo are well com- 
patible with this expectation. 

A fit with free /3+ gives a value of (3+ = 0.5(3) for 
amo = 0.04, and P+ = 0.5(4) for amo = 0.02. The 
probably overestimated error for AE± results in an over- 
estimated error of f3+. In fact our data do not have the 



necessary quality to investigate the scaling (4.11) with 
a free exponent Therefore, we fit them with fixed 
/3+ = 0.5. For j3pc{L,amo) obtained in this way we as- 
sume scaling as in the Higgs model EJ with vj^ = 0.5: 



Ppc{L, amo) - Pdamo) cx L ^^"^ 



(4.12) 



Our data are compatible with this assumption (Fig. g). 

The so determined points (3c{amo) of the EqoE line 
are listed in table || and plotted in Fig. |[ The point 



0.02 
0.04 
0.06 
0.10 
1.00 



0.654(7) 
0.671(4) 
0.684(4) 
0.707(3) 
0.80(2) 
0.848(4) 



0.304(5) 
0.296(4) 
0.292(3) 
0.285(3) 
0.28(1) 
0.263(2) 



TABLE II. Estimate for the position of the critical end- 
points Eamo on the EooE line (L+) in the infinite volume 
hmit. 



at amo — 1-00 was only estimated and the uncertainty 
included in error bars |lq ]. 

These results for Pd^mo) are extrapolated into the 
chiral limit. The curves should approach the symmetry 
plane with the critical exponent ^Aj: 

(3c{amo) - f3c{E) cx A,(amo) cx (amo)"^^* (4.13) 

A fit of the data obtained for amo < 0.10 gives (j)At « 1.6 
(Fig. |l^). This value is in agreement with that obtained 
by means of the relation (4.S) from the values of the 
exponents determined in the next section. There we find 
(f>At = 1.8(1). The extrapolated (3 value for the point E 
is /3e — 0.625. Of course, with three free parameters used 
to fit four data points the error is large and uncertain. 

The satisfactory quality of the fit and the agreement 
of both methods for the determination of (j)At indicates, 
that we may actually overestimate the errors in the whole 
procedure. E.g. fixing (jjAt — 1.8 in eq. ( 4.13| ) reduces 
the error for /3e without reducing the quality of the fit 
shown in Fig. |l^ and gives /3e — 0.62(1). 

In summary, we estimate the coordinates for the point 
E in infinite volume to be: 



/3e = 0.62(3) 



(4.14) 
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KE = 0.32(2) 



(4.15) 
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FIG. 9. Projection of the critical line of endpoints EooE 
and the line of triple-points TToo onto the /3-amo-plane. (In 
Fig. |l| this corresponds to a view from below for amo > 0.) 
Shown is also the critical line NE of chiral phase transitions. 
The points have been determined on G"* and 8* lattices. The 
error bars reflect the uncertainty of the extrapolation into the 
infinite volume. 



Ai= 1.6(13) 
E= 0.625(47) 
- A= 0.35(24) 




FIG. 10. The values of Pdarno) for amo < 0.10 ob taine d 
from extrapolation into infinite volume by means of (4.12). 



The fi t shows their extrapolation into the chiral limit using 



The rather small improvement of the precision compared 
to our earher pubhcation 0| shows how difficuh the de- 
termination of the position of the tricritical point is, if no 
simulations in the symmetry plane are possible. Never- 
theless, our present determination is much more reliable 
due to the use of the scaling analysis. 



V. CRITICAL AND TRICRITICAL EXPONENTS 



A. Exponents ft and /3„ 



iM3 



We have seen in section [II B that the scaling behavior 
at the point Eoo is mean-field-like. We now repeat the 
analysis by means of Fisher zeros also for small fixed 
aniQ and determine Vamo: the value of v at fixed amg, 
pursuing two aims. We want to check the universality 
along the EoqE line, and we want to look for a possibly 
different scaling behavior in the immediate vicinity of the 
tricritical point E. 

To use the method of the Fisher zeros, one coupling has 
to be fixed. For each aniQ we have fixed P at the position 
of ^amo ■ At this P we have done simulations for different 
K and afterwards determined the zeros in the complex 
K plane. The accuracy of the method is limited by the 
uncertain precision of (3c{amQ), as given in table ||. To 
estimate this dependence, we have made measurements 
at three /3-values for aniQ = 0.02. Between 9600 and 
64000 HMC-Trajectories have been done for each P value 
at 5-7 K-points. 

Fig. ^ shows the scaling behavior of the imaginary 
part of the Fisher zero in k for amg = 0.04 and 0.02. In 
all cases, no deviations from a linear behavior could be 
observed. The value 1)+ = i>o.o4 = 0.328(4) for aniQ — 
0.04 is in excellent agreement with the result at aniQ = 
oo. For amo = 0.02 it was more difficult to find the 
critical coupling /3c(0.02). From the three measurements 
we estimate ?'o.02 = 0.33(4). 

Investigation of the cumulants yields results for a^i- / 1>+ 
which are consistent with the use of the Josephson rela- 
tion. This indicates that the hyperscaling hypothesis is 
fulfilled. The calculation of i>+ with use of the Joseph- 
son relation yield consistent results with a little bit larger 
error. 

The finite size scaling theory above the critical dimen- 
sion should be applied with care. It is possible that in 
spite of consistent scaling it is not the exponent v which 
is observed^. In our case the values of iy^ obtained by 
this method and from the scaling of the fermion mass are 
consistent, however. 



■^We thank K. Binder for a discussion on this point 
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fit: A'L"""' 

*:,3 = 0.67 
5= 0,3282(37) 
A= 0,945(54) 
X?„=0,028390 
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am(, = 0.02 

fit: A'lr"': 

^:/3 = 0.64 
1i= 0,3749(88) 
A= 0.653(61) 
X/a.,= 4, 028770 

□ :^ = 0,655 
v= 0,3079(49) 
A= 1,32(11) 
X/a„,= 0, 000153 

*:,S = 0,66 
5= 0,2912(41) 
A= 1,69(15) 
X5„=0, 975449 



FIG. 11. Scaling behavior of the edge singularity: (a) for 
arriQ = 0.04 and (3 = 0.67, (b) for aniQ = 0.02 and 13 = 0.64, 
0.655, 0.66 on quadratic lattices. (3 = 0.67 is our current best 
estimate for the position of the endpoint Eamo at amo = 0.04. 
For arriQ = 0.02 we have investigated three different /3 values, 
[3 — 0.655 is nearest to the endpoint. 



Therefore, we interpret our results as a good confirma- 
tion of the universality along the EooE fine. The mea- 
sured values are nearly identical to those at amo — oo. 
Also the logarithmic corrections seem to tend into the 
same direction. Since we found similar results also for 
/?+ it is likely that all exponents along the EoqE line are 
independent of amg. 

Assuming that for the small values of amo, we could 
investigate, the dominance region of the tricritical point 
is already achieved, we expect that for amo ^ the 
exponent turns over into i^u — t^t- This implies that 
the subsidiary exponents with the index u are identical 
or at least very similar to the exponents along the EqoE 
line, which are the corresponding classical values. The 
measured value then corresponds to t't = i. This means, 
that the value is different from the classical values of a 
tricritical point. The corresponding classical prediction 

= ^ is hardly compatible with the data. 

A crossover to the exponents of the tricritical point 
at small amg could be expected also for the exponent 
/3+, which we measured along the EqoE line (in the S+ 
plane) in section IVC. However, also here we could not 
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FIG. 12. as function of amo at P — 0.55 on 6*, 
lO^'-lattices for different k. 



and 



observe any amo dependence and /3+ is compatible with 
the mean field exponent of a critical line /3+ = i down 
to amo = 0.02. We therefore interpreted this results as a 



indication, that also /?„ is close to 



2 ■ 



B. Estimate of St from R-^ 

As suggested in |^ it is possible to determine the 
'magnetic' exponent S for fermionic theories by measur- 
ing the susceptibility ratio R^^ for different small amo 
around the critical point. We have measured i?^ for dif- 
ferent K values. Inside the broken phase (k < Kc), we 
expect a curve which approaches i?^ = for amo 0. 
In the symmetric phase (k > Kc), we expect the curve 
to approach i?^ — 1 for amo ^0. At the critical point 
(k — Kc), the line should be horizontal for small amo and 
the corresponding value of Rtt is equal to 1/5 — l/St- 

We have tried this method in our model at /3 = 0.55 on 
the NE fine (Fig. |l^). For k = 0.37 the curves bend down- 
ward when approaching the chiral limit, for k = 0.38 
they bend upward. This indicates, that the critical k is 
between 0.37 and 0.38, in agreement with 0.376(5), our 
estimate based on the modified gap equation [Q. The 
estimated horizontal line separating both phases is at 
Rtt ~ 0.30(5). This is in good agreement with the mean 
field value of <5 = l/R^ = 3. Both results confirm our 
earlier result that the phase transtion a,t /3 — 0.55 is 
mean- field-like 0. 

Our results for (3 = 0.64, close to the tricritical point, 
are shown in Fig. |l^. A horizontal curve is found for 
Kc = 0.307(2). It extrapolates to R.:,{kc) = 0.10(5). We 
estimate that St = 10^4° under the assumption, that 
the increase of the volume does not change this picture 
significantly. Indeed, the data for the 8''' and 10^ differ 
only very little. 

To estimate the sensitivity of these results on /3 we 
have obtained lower statistics data also a,t f3 — 0.62, our 



12 



VI. SPECTRUM 




FIG. 13. Rtt as function of amo close to the tricritical point 
at /3 = 0.64 on 8^ and 10^-lattices for different k. 



current best estimate for /3e- The results are very similar 
to those at (3 — 0.64. This indicates that for (3 = 0.64 we 
are already in the influence region of the tricritical point. 

Also an explorative investigation with the method of 
the Lee- Yang zeros in the complex amo plane, as we 
applied recently in the three-dimensional model p6| , p^ , 
confirms the large value for St. Using the lattices 4^ 
and 6'^ we get for (3 — 0.64 and k — 0.307 the estimate 
D ~ 0.275, which corresponds to St — 10. 



To analyze the possible physical content of the contin- 
uum limit taken at the point E, it is important to study 
the spectrum. This investigation has the advantage that 
it is relatively independent of the missing theory of tri- 
critical points in four dimensions. Here we give only the 
most important results. A tabular overview of the mea- 
sured values (which could be presented even graphically 
only in part) can be obtained from the authors. 

Most of the shown results have been obtained for fixed 
K = 0.30. The reason for this is the observation, that 
the shift of the endpoint Eamo with aniQ is smaller in 
the K direction than in the (3 direction. The advantage 
of this is, that for different amo the endpoint Eamo can 
be approximately hit with one k. The difference of this 
chosen value of k from our current best estimate ke = 
0.32(2) is due to the underestimated remaining shift at 
the beginning of the large scale simulation. The value 
K = 0.30 corresponds the the position of the endpoint 
Eamo at amo ~ 0.03. 

Some measurements with less statistics have been per- 
formed for /? = 0.64 and k = 0.305. They confirm the pic- 
ture presented here. In particular, they indicate a com- 
mon scaling behavior in a whole region around E, e.g. 
independent of the direction of approach to this point, 
provided it is not tangential to the critical lines. 



C. Summary of results for the tricritical exponents 

We have determined three independent exponents at 
the tricritical point: 



f3u 
St 



0.33(5), 
0.5(2), 

10^ 



-10 
-4 ■ 



(5.1) 
(5.2) 
(5.3) 



These values disagree with the classical values for tricrit- 
ical points expected in four dimensions, i^t — 0.5, /?„ = 1, 
St = 5. 

The errors take into account only uncertainties in the 
measurement. We cannot estimate possible systematic 
errors. In particular we have assumed that at amo — 0.02 
on the used lattices we observe the asymptotic scaling 
behavior of the tricritical point. Some support for this 
assumption is obtained in the spectrum analysis in the 
next section. It is plausible also because typically tri- 
critical points have a large region of influence and the 
corresponding deviations from the scaling behavior are 
small. 

If one assumes the validity of the scaling laws, all fur- 
ther exponents are fixed. We only could check with good 
precision the Josephson relation between at and Vf Some 
very crude check was possible for u^, and another one for 
Pt is described in Sec. VI A4. We found no indications 



A. Fermion, boson and gauge-ball spectrum 

1. Uniform behavior as a function of amp 

The pictures showing the behavior of the mass of 
fermion and scalar boson S in the vicinity of E have been 



shown already in section IV B. The very similar position 



of the pseudocritical area of ams and amp for different 
volumes and small amo for k = 0.30 ~ ke suggests to 
express one mass as a function of the other. 

As amp is monotonously decreasing with increasing k 
or f3, and is well measurable, we plot ams as a function 
of am p (fig. ^4|) . In this figure the broken phase is to the 
right and the symmetric phase to the left. The border 
between both phases is at the minimum of the boson 
mass. In the infinite volume and the chiral limit, the 
symmetric phase would reduce in this plot to one point, 
amp = 0. This kind of plot as a function of amp is 
possible for all couplings near to the NE line, but ams 
scales only at the point E. 

To extract continuum physics in this way, it is neces- 
sary to check that all curves are close to a uniform curve, 
when the tricritical point is approached. In principle, 
this means the fourfold limit ^ — *■ oo, amo — *■ 0, k — > ke 
and P Pe- 



for a violation of the scaling laws. 
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^ = 0.30 




FIG. 14. Boson mass ams as function of the fermion mass 
amp- The broken phase is to the right from the minima. 
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FIG. 15. Mass ratio ms/niF as function of the fermion 
mass. 



One can see that in the broken phase, ams ^ a func- 
tion of amp is nearly independent of volume and am^. 
For other k ~ k-^ and P ~ (3-^ nearly the same curves re- 
sult. From this nontrivial observation we conclude, that 
our data for fixed amg and n = 0.30 correspond approxi- 
mately to a path towards E for amo = 0. For small amp, 
at the transition into the symmetric phase, the boson 
mass ams increases rapidly. As expected, the minimum 
shifts to the left for increasing lattice size and decreasing 
amo, which corresponds to the approach to the critical 
theory (chiral limit in the infinite volume). 

In the following, all observables are plotted as a func- 
tion of amp. To make the figures more clear, usually 
only the data on the 8'^24 lattice are shown, as long as 
the effects of finite volume in this kind of plot are small. 



2. Scaling behavior of ams 

To investigate the scaling behavior of ams we look 
at the mass ratio ms/mp. As shown in fig. |l5|, this 
ratio is nearly constant, and close to 0.5 in the broken 
phase. A small decrease of this ratio for decreasing amp 
is indicated. Most probably the reason for this is, that for 
nonzero amo, the boson mass ams vanishes at the point 
Eomo in the V ^ oo limit, whereas amp stays finite. We 
have checked that approximately the same mass ratio 
is obtained also on other paths into the point E. This 
strongly suggests that this mass ratios is preserved also 
in the continuum limit. This would mean that the scalar 
boson would survive with approximately half the fermion 
mass. 

This similar scaling behavior means, that both observ- 
ables have a common critical point. This is our strongest 
argument for the coincidence of chiral and Higgs phase 
transition at the point E which is thus a tricritical point. 

The fact that this nice scaling behavior can be ob- 
served for all investigated amo < 0.04 further suggests 



that these values of the bare mass are inside the dom- 
inance region of the tricritical point. In the influence 
region of the EooE critical line amp would stay finite 
whereas ams would scale to zero. This should at least 
show up in scaling deviations. 



3. Gauge balls and vector boson 

We also measured the mass of two gauge balls. 
The scalar gauge ball with the quantum numbers 0++ 
(Fig. |l^ has nearly exactly the same values as the scalar 
boson S obtained from the cf)^ — (j> correlation function 
(^.16| ). We have also measured the cross-correlation be- 
tween both channels at some points and found a mass 
in good agreement with that from individual channels. 
This strongly suggests that in both channels we see one 
state, which can be interpreted both as a scalar boson 
and gauge ball. In the Nambu phase, the second inter- 
pretation might be more natural as it holds in the whole 
phase, including k = 0. Nevertheless, we continue to 
denote the state by S. 

We also looked at the gauge-ball channel with the 
quantum numbers 1"' , but we could not observe any 
light particle near the tricritical point. 

The mass of the vector boson amy does not scale 
(Fig. |l^. The mass decreases significantly in the sym- 
metric phase but stays large at the phase transition 
(amy > 1). 



4- Scaling of amp and (xx) 

We have tried to estimate the tricritical exponents also 
from the fermion mass amp and the chiral condensate 
(xx)- The proper method would be to extrapolate the 
data first to the infinite volume and then into the chiral 



limit. As discussed in section IVB, we are not able to do 
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FIG. 16. mass ration mc/m-F of the scalar gauge ball 
(0^"'') as function of the fermion mass. 
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FIG. 17. Mass of the vector boson amy as function of (a) 
/3 and (b) amp for k — 0.30 on the 8'^24 lattice. The phase 
transition (minimum of ams) is marked in (a) with arrows. 



FIG. 18. Scaling behavior of amp (filled symbols) and (xx) 
(open symbols) as function of the coupling (5 for k = 0.30 on 
8^24 lattices. For each amo a simultaneous fit to (xx) ''■iid 
amp with one common critical coupling was done. The fit 
interval is marked with arrows. The results of the fit are 
given in the legend (the amplitudes are not shown). 



(as determined by the minimum of ams) in the broken 
phase. At these points the value of the condensate is 
still around or above 0.3. So we end up with 6 measured 
values and 5 free fit parameters. 

The so estimated parameters (3pc, Pt and vt are listed 
in Fig. |lj. The values of f3pc are in nice agreement with 
the minimum of the boson mass ams- The results for 
critical exponents are in rough agreement with the values 
of i^t = 0.33(5) and Pt — 0.12tg^ obtained in section 
0. (T he latt er va lue is obtained from (x2) by means 
of ( |4.9| ) and ( 4.10| ).) These results support the values 
of the exponents (eq. ( ^.l| ) and (5.3)) obtained by other 
methods at the tricritical point. 



this. Therefore we have investigated whether both effects 
can be absorbed into the aniQ and volume dependence of 
the pseudocritical coupling (3pc. This is suggested by the 
fact that the value of ams /amp is nearly independent 
of amp and of the used (small) am^. As ams scales 
for amo ^ 0, this should be approximately so also for 
amp and (xx)- (^"^s itself is less suitable for such an 
investigation, because of the larger errors). 

Correspondingly, we plot the results for the fermion 
mass amp and the chiral condensate (xx) for different 
amo as a function of the coupling /? (Fig. |l^) . We expect 
the approximate scaling behavior 

amp = Ap{(3p, - (6.1) 
{xx)^MPpc-Pf' , (6.2) 

where Ppc depends both on the volume and amQ. To get 
stable results we did, for each amo, a fit with a common 
pseudocritical couphng Ppc- We choose for the fit of both 
observables the 3 points closest to the critical coupling 



B. Properties of tt- meson 

The mass of the tt- meson (0 ^) can be measured very 
reliably. In the broken phase we expect the validity of 
the PCAC relation 

(am^r)^ oc amo • (6-3) 

As shown in Fig. |l^, this relation is nearly perfectly ful- 
filled for P = 0.64, K = 0.305 (near to the tricritical point 
in the broken phase). Thus it has the expected scaling 
behavior of a Goldstone boson. For k = 0.31 (transi- 
tion to the symmetric phase) an expected deviation is 
observed. 

To investigate the fiavor symmetry restoration we com- 
pare m Fig. 1^ the mass in the channels (1) and (2) which 
belongs to the quantum numbers ^ (table I in [Q). 
These masses are labeled am^*\ with the corresponding 
channel in the exponent. In channel (2) also the first 
excited state could be measured and is labeled am^ 2- 
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FIG. 19. (am,r)^ from channel (2) as function of amo. The 
PCAC relation ( p.3| ) corresponds to a straight line through 
the origin. 
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FIG. 21. Pion decay constant a/^ as function of amp for 
K = 0.30. 
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FIG. 20. Mass of am^ from channel (1) (circles) together 
with mass am,r and first excited state of channel (2) (squares 
and triangles, respectively) as function of amp- In channel (1) 
a simultaneous fit to the data with point source and smeared 
source has been used. 



A light particle can be observed only in channel (2). 
This is the Goldstone boson, which scales corresponding 
to the PCAC relation. Its mass is different from that in 
the channel (1). The first excited state in channel (2) 
is close to the mass in channel (1). Both scale with ap- 
proximately twice the fermion mass. The small observed 
deviations for small amp are in the symmetric phase. 

We cannot see restoration of flavor symmetry with par- 
ticipation of the massless Goldstone boson. But an agree- 
ment of the massive contributions with the corresponding 
quantum numbers seems to show up. This behavior of 
the Goldstone boson, which is massless in the chiral limit, 
and its role in the flavor symmetry restoration is not yet 
understood. Similar behavior was observed e. g. in the 
NJL model (21 1. 

We note that the data indicate presence of a light pseu- 
doscalar in the symmetric phase, as seen e.g. in Fig. |2^ 




- am„ = 0.02 



(3 = 0.00 (6=16) 
,3 = 0.55 (8=24) 
,« = 0.84 (8='24) 
^ = 0.30 (8=24) 



FIG. 22. Pion decay constant a/^ as function of amp for 
amo ~ 0.02 and different (3 and k, respectively, along the NE 
line. 



at small amp . In the chiral limit the massless Goldstone 
particle seems to change at the phase transition into a 
(bound?) state of two massless fermions F. 



1. Pion decay constant 



Fig. ^ shows the pion decay constant a/^ (eq. ( ^.10| )) 
as a function of am p . The value of a/^ is nearly indepen- 
dent of amo s-iid for larger amp also nearly independent 
of the volume. For small amp there is a clear tendency 
for decrease of af.^ on larger lattices. 

Fig. shows a comparison of af-^ for different f3 along 
the NE line. The smallest values for a/^ can be observed 
close to the point E, but the difference between (3 = 0.55 
and (3 = 0.64 is very small. 

If the pion decay constant o/tt scaled like a mass, 
its plots as a function of amp should give a strait line 
through the origin. This is not observed, although on 
larger lattices there is a slight reduction observable even 
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FIG. 23. p meson mass from (a) channel (4) and (b) chan- 
nel (3) for K = 0.30 on the 8^24 lattice. The line shows the 
2amF threshold. 




amp 

FIG. 24. a meson mass for channel (1) for k = 0.30 on 
the 8'^24 lattice. The full line shows the 2amF threshold, the 
three dotted lines twice the tt meson mass. 



in the broken phase. We do not know, if this is a hint 
towards a larger finite size dependence. One may expect 
such a large dependence on the basis of estimates made 
by means of the Schwinger-Dyson equations for the NJL 
model (Fig. 33 in Ref. ||2^]). At present we cannot ex- 
clude that /tt/wf diverges in the continuum limit. This 
could be an indication of the trivial continuum limit [0 . 
In any case the largest value fTr/mp ~ 1/3, we found in 
the broken phase around amp ~ 0.4, is a lower bound 
for this ratio. 

C. Further mesons 

1. p meson mass 

The mass of the p meson (1 ) can be measured quite 
well. We use a fit to the propagator with a smeared 
source and sink, which suppresses excited states quite 
well. We get the same mass also if we use the point source 
and sink and do a fit with two states with negative (and 
one with positive) parity. 

The interesting feature of the p particle is the fact, that 
we can observe it in two channels. As is shown in fig. 
the results are nearly identical, especially near to the 
phase transition {amp ~ 0.5). So the flavor symmetry 
seems to be restored within good precision, at least for 
the p. 

The values close to the phase transition are nearly ex- 
actly equal to 2am p (line)|^ and the p meson thus scales 
in excellent agreement with the fermion F. We cannot 



^The somewhat larger values for arup in [|l^ are due to an 
insufficient consideration of the excited states. 



distinguish, whether the p meson is a bound state or a 
resonance. 

2. a meson mass 

Fig. [2^ shows our measurement of the <j meson (0++) 
mass without consideration of the annihilation part. This 
simplification is necessary due to computer time restric- 
tions. It is striking, that the a meson mass is especially 
for small otoq nearly independent of amp. In the broken 
phase it does not increase with it. We also observe some 
finite size effects, with tendency of increasing am^ with 
increasing volume. Furthermore, as seen in Fig. the 
dependence on the bare mass is large. The data therefore 
do not allow us to extrapolate m^ to the infinite volume 
and chiral limits. 

It is probable that the a meson decays in two tt mesons. 
Therefore for each amo the dotted curve shows 2am7r. 
All measurements are within the error bars on or above 
this threshold. Possibly, what we observe is not the a 
resonance but two tt mesons. It should be noted that 
also in the NJL model the cr meson shows large finite 
size effects and dependence on the bare mass 

The measurement of the mass of the a meson would 
have been very interesting, as it is the most natural candi- 
date for the Higgs boson. Because also the scalar boson 
S has the same quantum numbers, we would have ex- 
pected a mixing of those two states. This does not seem 
to be so. Thus we have to leave open the question, with 
which of the scalar particles observed on the lattice the 
(composite) Higgs boson has to be identified and what 
mass it has. 
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VII. SUMMARY AND CONCLUSION 

The tricritical point E in the xU(t> model turns out 
to be a very complex phenomenon in four-dimensional 
quantum field theory, presenting numerous challenges. 
No reliable way is known to study it analytically. Its 
numerical investigation faces tremendous obstacles: the 
analysis of numerical data is made without a plausible 
analytic scenario, two couplings have to be fine-tuned, 
the use of dynamical fermions is necessary, and the chiral 
symmetry plane is not yet accessible to simulation. The 
last two obstacles make the investigation of this point 
substantially more difficult than the study of tricritical 
points in metamagnets and other systems in statistical 
mechanics. Though statistical mechanics provides the 
conceptual understanding of tricritical points which we 
have heavily used, there is actually no experience with 
tricritical points in four dimensions to compare our re- 
sults with. 

For these reasons our conclusions can be only tenta- 
tive. Nevertheless, they justify the effort^ we have spent: 
there is fair chance that the tricritical point E exists and 
defines an interesting quantum field theory in four di- 
mensions. Though it has various features analogous to 
the NJL model, the gauge field plays an important dy- 
namical role. It may be very complex and completely 
unaccessible by perturbation theory. But it may provide 
a new alternative for dynamical explanation of the masses 
of fundamental constituents in the standard model. 

The simplest way to justify this hope is to describe 
our results as a microscopic model for a strong Yukawa 
coupling, postulated in the standard model e.g. for the 
top quark. We have done this in a separate letter |l^] 
^. In the vicinity of E, the Yukawa coupling between F 
and TT emerges naturally as the Van der Waals remnant 
of the very strong interactions of the xU4' model. 

The most surprising result is the clear signal for the 
scalar S of mass ms — l/2mF- It can be interpreted ei- 
ther as a bound state of the pair of fundamental scalars 
4>^,(j), or as a gauge ball. Both channels mix strongly. 
On the other hand, the would-be Higgs boson, the fun- 
damental fermion-antifermion bound state a = xx, is 
elusive. It does not show scaling properties allowing an 
extrapolation to the continuum limit. It does not mix 
appreciably with S. 

The most perplexing result is the value of the tricritical 
exponent Vt — 1/3, which is a nonclassical value. Could 
it mean that the continuum theory is nontrivial? Stan- 
dard lore in statistical mechanics is that in and above 



*E.g. equivalent of 2 years on a CPU of one Gflop speed, 
producing about 10* HMC trajectories on various lattices. 

^Some important technical details are given in the Appendix 
of this paper. 



three dimensions tricritical points are classical [g| . How- 
ever, this is based on the experience with spin systems 
and scalar fields. Strongly coupled gauge theory with 
fermions might be different. Thus point E is a challenge 
also for statistical mechanics. 

We feel that our aim to understand the tricritical point 
in the xUcj) model was a little bit ahead of time. Though 
we used the most advanced methods, they were not pow- 
erful enough. The computational resources should be 
also larger by 1-2 orders of magnitude. Algorithms for 
simulation in the chiral symmetry limit are needed. Last 
but not least, the interest in a search of strongly cou- 
pled theories beyond the standard model should be much 
higher than the current widerspread beliefs. 
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APPENDIX A: DETAILS ABOUT THE USED 
OPERATORS 

1. Smeared sources for meson propagators 

To reduce the contribution of excited states in the me- 
son propagators we have implemented gauge invariant 
smeared sources |Q. We are not aware of any such an 
implementation for staggered fermions. It requires that 
the even-odd separation is preserved. Therefore we have 
transported the source with two link term to the next to 
nearest neighbour. This reads as 

1 — Da 

< (l){x, t) + a ^ [Ui{x, t)Ui{x + el, t)4>{x + 2el, t) 

-hf/*(x- el,t)t/;(f - 2el,t)</)(f - 2el,t)] |. (Al) 

We have chosen values a = 0.01 and a — 0.02 and 20 
smearing iterations. The resulting smearing radius was 
R ~ 0.89 and R ~ 1.47, respectively. The latter gave the 
better results. Compared to QCD these radii may seem 
very small, which might be due to larger masses in our 
case, however. 
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2. Effective Yukawa coupling yn. 

We have done measurements of yj^ in the momentum 
space, similar to 27 1. The meson- fermion three point 
function is 



q')^^{j) = i_ g-iP4tl +194*2 g-ipXj+ilJ£2 

*2 -£11312 

M„(p-q,y4)0^(ai2,i2)x(£2,i2) ) , (A2) 



with 



Ma(po, y4) = 5] e'P^°V,(y)x(y, y4)x(y, 2/4) , (A3) 



and <fa{y) = 1 for the a meson (scalar), and '^■n{y) = 
Viiif) = (— I)!'i+y2+y3 for the tt meson (pseudo-scalar). 
The underscore of Xi and X2 indicates that the spatial 
part of the vector has even coordinates. The correspond- 
ing sums run over the 2^ cubes of the spatial lattice. The 
measurements were done for the momenta 



Pi = 



54 for the a meson 
q4 + t: for the n meson 



94 

2n 



(p,g) = (0,0), (_el,0),(-el,^el). 



±y (A4) 
(A5) 



This choice of p4 guarantees, that only states of the right 
parity contribute. For 54 only the smallest possible values 
are considered and the results are averaged over both 
signs of q^. The different spatial momenta are evaluated 
separately. The spatial momentum of the meson is po = 

p- q- 

For the implementation we neglect the disconnected 
parts and write by means of the fermion matrix My^ 



Gi"\p,q) 



1 

T 



^e^(P-^)Va(y)e(2/) 



^M-^\0(xi)e- 



-ipxi 



^ M^-40(x2)e-(^+(°>))- \ ), (A6) 



where the four- vectors xi = {xi,ti), y — {y, y^), ... have 
been introduced. In this notation it is well noticeable, 
that for our 6 different momentum pairs {p, q) for both G3 
all together 8 matrix inversions are needed. Due to the 
source on the whole lattice the signal is very good and we 
have done a measurement only on each 8'^ configuration. 



The effective Yukawa coupling y^-* is now obtained 
from the comparison of the Monte-Carlo data for the 
three-point function, with that obtained in the tree level 
approximation of an effective lattice action, which de- 
scribes the interaction of the staggered fermion fields F, 
F and a (pseudo-)scalar field <I> with a coupling term of 
the form 



Vr Y1 '^ix)F{x)F{x) (scalar) , 



(A7) 



(TV) 



'^e{x)^{x)F{x)F{x) (pseudoscalar). (A8) 



The connected part of the three-point function is 
,(") 
T 



Gt^ (p, 9) - - ^ ^ E ('^)^" - ^) 



J E Vii^y'e<P+<^°^''^>'+^^GFixi + u},0)A (A9) 



with 



1. 



(AlO) 



Here uj runs over corners of the elementary three- 
dimensional cube. In the tree-level approximation, Gq, 
and Gf are the free propagators for the meson and the 
fermion, the tilde indicating the Fourier transformation 
of the meson propagator. These propagators are replaced 
by the full propagators from the simulation and the wave- 
function renormalization constants Zf and are in- 
cluded: 



GaiPo) = V Za{po)Oa ^ ' 



^(x)x(x)x(x) ^(0)x(0)x(0)) (AH) 

(A12) 



GFix,0)p-= ^ZF{p){xix)xm , 



with 9cr = —I and 6'7r = 1 to correct for the negative sign 
of the a propagator. 

After these replacements and identifications of 



G3 ip,q) with the measured values of expression (A6) 



we obtain y^\p,q) from eq. (A9). j/r, should be real 
and only slightly dependent on the momenta. For the 
determination of Z{p) we have measured the correspond- 
ing propagators at momentum p and then performed a 
fit with the free propagator. 

Using these definitions we have measured the real and 
imaginary part of j/r . The implemented momenta com- 
binations have been numbered corresponding to equation 
(^) from 1 to 3. 
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FIG. 25. Results for the real and imaginary part of 
3/R of the TT meson for different momenta at k = 0.30 and 
0.02 on the 8^24 lattice as a function of amp- 



For the effective coupling of the tt meson we get a con- 
sistent picture (Fig. p5|) . The imaginary part of t/R is very 
small, and for different momenta we get approximately 
the same values. Because of this agreement we restrict 
ourselves to the evaluation of the data with vanishing 
momentum (No. 1). 

For the a meson we failed to get a reliable j/r. The 
imaginary part is not really small and the real part is 
only for the momentum combinations 1 and 2 (at least 
one fermion momentum vanishes) approximately equal. 
The inconsistencies are larger in the broken phase. The 
problems might be related to the neglected disconnected 
parts but also to the fact, that the a meson is probably 
only a broad resonance. In this sense this measurement 
shows again the problems we met already in the mea- 
surement of the (7 meson mass. 
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